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ABSTRACT 

The acoustic radiation from isotropic turbulence is computed numerically. A hybrid 
direct numerical simulation approach which combines direct numerical simulation (DNS) 
of the turbulent flow with the Lighthill acoustic analogy is utilized. It is demonstrated 
that the hybrid DNS method is a feasible approach to the computation of sound generated 
by turbulent flows. The acoustic efficiency in the simulation of isotropic turbulence 
appears to be substantially less than that in subsonic jet experiments. The dominant 
frequency of the computed acoustic pressure is found to be somewhat larger than the 
dominant frequency of the energy-containing scales of motion. The acoustic power in the 
simulations is proportional to eMf where e is the turbulent dissipation rate and M t is the 
turbulent Mach number. This is in agreement with the analytical result of Proudman 
(1952), but the constant of proportionality is smaller than the analytical result. Two 
different methods of computing the acoustic power from the DNS data bases yielded 
consistent results. 
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in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 24681. 
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Introduction 

In 1952, Lighthill posed the problem of estimating the sound radiated by “a fluctuating 
fluid flow occupying a limited part of a very large volume of fluid of which the remainder is 
at rest” - the fluid flow of course is either ordered, with coherent distributions of vorticity, 
or disordered and turbulent. It is only now, four decades later, that the simplest possible 
model of this problem i.e., sound radiated from a limited body of isotropic turbulence 
bounded by a quiescent fluid, is amenable to numerical simulation. 

Isotropic turbulence in a bounded domain is a model wherein the turbulence is unaf- 
fected by the boundaries enclosing the fluid, and furthermore the statistical moments are 
spatially invariant and independent of orientation. Isotropic grid turbulence is a similar 
idealization, in that the turbulence is enclosed by wind tunnel walls and the homogeneity 
of the turbulence in the central region is known to be unaffected by the wall boundary 
layers. The problem of sound emitted by isotropic turbulence has been investigated by 
Proudman (1952), and more recently by Lilley (1993) in an analytical study and by 
Witkowska and Juve (1993) in a large eddy simulation (LES). The isotropic turbulence 
is represented in our numerical study by a periodic box which contains the fluid motion 
and is surrounded by a layer in which the velocities decrease to zero. The fluid motion 
inside the periodic box is simulated, and it is assumed that the sound radiated by the 
periodic box is unaffected by the surrounding thin layer. 

Two possible numerical approaches to the prediction of sound produced by turbulence 
are the direct numerical simulation (DNS) and the hybrid DNS. In a DNS, all relevant 
scales of motion are numerically resolved in a solution of the compressible Navier-Stokes 
equations; while in a hybrid DNS only the turbulent flow is resolved numerically. Sarkar 
and Hussaini (1993) have estimated the turbulence Reynolds numbers Re t that could 
be obtained for a three-dimensional calculation of sound generated by turbulence with 
the two methods for a given spatial grid. It was found that for Mt = 0.01 and 0.1, the 
Re t that can be simulated by DNS of sound on a uniform grid is respectively 0.004 and 
0.08 times smaller than that achieved in the hybrid DNS method. Evidently, the much 
larger computational volume required in the DNS of sound approach leads to a much 
lower turbulence Reynolds number relative to the hybrid approach. By stretching the 
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computational grid or by using overlapping grids, the Re t obtainable with the DNS of 
sound approach can be increased. However, the effect of grid stretching on the accuracy 
of the far-field acoustic pressure which is several orders of magnitude smaller than the 
aerodynamic pressure has to be carefully examined. Crighton (1993) points out that 
computations of low Mach number aeroacoustics by the DNS of sound approach have 
to contend with a number of difficulties which include the large disparity between the 
acoustic length scale and the turbulence integral length scale, the much smaller energy in 
the acoustic far field relative to the turbulent energy, and the multipole structure of the 
basic acoustic source. Due to these reasons, Lighthill (1993) recommends the use of an 
acoustic analogy for the computation of sound radiated from low Mach number flows. The 
present work employs the the hybrid DNS approach with the Lighthill acoustic analogy 
for the purpose of computing the sound generated by three-dimensional, fully turbulent 
flow. Lighthill’s analogy has been used previously to compute the sound numerically from 
the large-scale instability waves in a shear layer by Gatski (1979), but we are not aware 
of any previous attempt to use the analogy within the framework of three-dimensional 
DNS of turbulence. 

The acoustic analogy which was first proposed by Lighthill (1952) reduces the aeroa- 
coustics problem to an inhomogeneous wave equation. In his analogy, Lighthill chose the 
homogeneous part to be a linear wave equation for the density and obtained the acoustic 
far-field as an integral of the equivalent acoustic sources in a uniform medium at rest 
which replace the fluid motion inside a bounded volume. Furthermore, by considering 
the statistics of the acoustic sources, Lighthill deduced the eighth-power law at low Mach 
numbers for the net acoustic power radiated from a turbulent region. 

Lighthill’s analogy was followed by other analogies (Phillips (1960), Powell (1964), 
Ribner (1964), Lilley (1974), and Howe (1975)). In the analogies of Phillips, Lilley and 
Howe, the homogeneous part of the equation for the acoustic variable is a nonlinear 
convected wave equation rather than the linear wave equation in Lighthill’s analogy, 
and as discussed by Goldstein (1974), Crighton (1975) and Ffowcs Williams (1977) has 
the advantage of explicit consideration of the convection and refraction of sound by 
turbulence. Powell and Howe rewrite the term u • Vu as x u + Vu 2 /2 in their analogies 
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which leads to an explicit vorticity-dependent term in the inhomogeneous part of the wave 
equation. The explicit appearance of the vorticity in the forcing term is convenient in 
unsteady flows such as the Von Karman vortex shedding in the laminar wake of a cylinder, 
but perhaps not so important in high Reynolds number turbulent flows where the flow is 
more incoherent and cannot be idealized as a simple collection of idealized vortices. It is 
also more convenient in the context of vortex methods (Leonard (1985)). Applications of 
the Lighthill theory are consistent with experimental data on noise radiated from subsonic 
jets running at ambient temperatures. However in supersonic jets (Seiner(1992)), it 
may be necessary to account for additional physical features explicitly such as Mach 
wave radiation, jet screech, and other shock-associated noise. In our work on noise 
generated from isotropic turbulence at low Mach numbers, Lighthill ’s acoustic analogy 
is the preferred representation and is the one followed herein. 


Computation of the Lighthill acoustic analogy 


The acoustic analogy of Lighthill can be written in the space - derivative form 



r/cp) 

r 


d y, 


( 1 ) 


but for a far-field location of the observer point with |r| >> |t/|, Lighthill showed that 
the integrand in (1) can be approximated by the second time derivative taken at the 
retarded time yielding 


Px ( x,l) = i ^^/f jj (y,i- r /c 0 )<iy, (2) 

where pa is the fluctuation of pressure relative to the ambient, Co is the ambient speed 
of sound, and r = |x — y|. In Eq. (1), density fluctuations are neglected in the source 
term and the approximation, T t jjj ~ d 2 {pu^Uj)j dyidyj, is utilized since only low Mach 
numbers are considered in the present work. Here, p is the mean density which is of 
course a constant equal to the ambient density of the surrounding medium which is at 
rest. In Eq. (2), T^y, t-r/co) denotes d 2 (pu i u :j )/dt 2 evaluated at retarded time t,-r/c 0 . 

The major advantage of the time-derivative form of the acoustic analogy over the 
space-derivative form is that the quadrupole nature of the sound source is directly evident 
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in Eq. (2) but not in Eq. (1). However, the computer storage requirements of the space- 
derivative form is smaller by a factor of 6 relative to those of the time-derivative form 
which could present a significant advantage in large-scale computations where computer 
memory is a constraint. In Sarkar and Hussaini (1993), the relative advantages of Eqs. (1) 
and (2) were investigated for a test problem with an imposed quadrupole. The retarded 
time effect was accounted for by the time accumulation method wherein the observer 
time t was approximated by t ~ [(r + cor)/At]At where [.] denotes the integer value 
function, and At denotes the time step used for time advancement of the flow. It was 
found that the number of time points required for the spatial derivative form is too large; 
turbulence of the quadrupole source type requires 0(1 /M?) time points per oscillation 
of the source. Due to the stringency of this requirement, we choose the time-derivative 
form, Eq. (2), over the space- derivative form, Eq. (1). 

Furthermore, Eq. (1) is not quite appropriate here because of the periodic boundary 
conditions used for the turbulence simulation. Because of the non-zero velocity at the 
fluid boundary of the periodic box, the dipole contribution to the sound from the surface 
of the flow domain is non-zero and would be included in a calculation using Eq. (1). 
The surface contribution is an artifact of the boundary conditions and is not of interest 
to us. Eq. (2), on the other hand, gives directly the quadrupole sound characteristic 
of a finite volume of turbulent flow embedded in a medium at rest. We ensure that the 
periodic boundary conditions do not unrealistically affect both the flow field and acoustic 
source by keeping the length of the computational domain much larger than the spatial 
correlation lengths of the velocity field u % and the acoustic source T,j. 

In the case of unforced isotropic turbulence, the acoustic source term, f Tij( y, t — r/co) dy, 
decays in time. Consequently, the acoustic pressure is statistically unsteady and time av- 
eraging is inappropriate for statistical analysis. Short-time averaging was used in Sarkar 
and Hussaini (1993) to compute the statistics of the acoustic pressure. But such a method 
has since proven to be unsatisfactory due to the time scale of the turbulence decay being 
of the same order as the time scale of the acoustic fluctuations. The results presented 
here are substantially different from the preliminary work of Sarkar and Hussaini (1993). 

In the present work, ensemble averages are used to compute the acoustic intensity, power 
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and frequency spectrum. 20 cases were run with the same parameters and the same 
initial statistics, but with different instantaneous initial conditions. For each case, the 
observation points were distributed on two concentric spheres in the far-field, 10 on each 
sphere. Due to the isotropy of the turbulence, the 20 different simulations with 10 dif- 
ferent observer points at the same distance from the center of the flow domain lead to 
200 different samples of the acoustic pressure which is a sufficiently large number to ob- 
tain acceptable statistics. The scaled acoustic power is obtained by averaging over 400 
samples. 

Flow simulation method 

The turbulent flow inside a cubical domain is obtained by solving the compressible 
Navier-Stokes equations by a numerical algorithm which was originally developed to in- 
vestigate compressibility effects in isotropic turbulence (Erlebacher, Hussaini, Kreiss, and 
Sarkar (1990), Sarkar, Erlebacher, Hussaini, and Kreiss (1991a)), and homogeneous shear 
turbulence (Sarkar, Erlebacher, and Hussaini(1991b), Sarkar (1992)). The homogeneous 
turbulence problem permits periodic boundary conditions in all three coordinate direc- 
tions. It is of course necessary that the length of the computational domain be much 
larger than the integral length scale of the turbulence for obtaining realistic flow fields. 
Spectral accuracy is obtained by using a Fourier collocation method for the spatial dis- 
cretization of the governing equations. A third order, low storage Runge-Kutta scheme 
is used for advancing the solution in time. 

Initial conditions have to be prescribed for «/, p, p and T. The initial velocity field 
is split into two independent components, that is, u/ = uj' + uf \ where each component 
has a zero average. The solenoidal velocity field uj' which satisfies V-u^ = 0 is chosen 
to be a random Gaussian field with the power spectrum 

E(k) = k 4 exp(— 2Ar 2 /A£j) (3) 

where k m denotes the wave number corresponding to the peak of the power spectrum. 
The compressible velocity uf which satisfies V x u c/ = 0 is also chosen to be a random 
Gaussian field satisfying the same power spectrum, Eq. (3). The power spectra of the 
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two velocity components are scaled so as to obtain a prescribed Unns = yujuj, an d 
a prescribed x = u Sns/ u rms which is the compressible fraction of kinetic energy. The 
pressure p 1 ' associated with the incompressible velocity is evaluated from the Poisson 
equation 


r~7 2 V 

V p = — pu ■ -u 


hj ' 


V 

3 ,' ' 


( 4 ) 


It remains to specify the initial values of the thermodynamic variables. The mean density 
p is chosen equal to unity, and p is chosen so as to obtain a prescribed Mach number 
Umis/ySr / P characterizing the turbulence. The fluctuating density p' and compressible 
pressure p c ' are chosen as random fields with the power spectrum determined by Eq. (3) 
and prescribed p and p^ s . The pressure then becomes p — p + p 1 ' p c> , the density 
is p = p + //, and the instantaneous temperature T is obtained from the equation of 
state p = pRT. We note that the splitting of velocity and pressure into incompressible 
and compressible components is confined to the specification of initial conditions. No 
such split is carried out while computing the acoustic sources during the time evolu- 
tion. Although the numerical algorithm is capable of simulating compressible turbulence 
with non-zero density fluctuations and non-zero dilatation, we consider low Mach num- 
ber, quasi-incompressible turbulence in the present paper. The acoustic source T,j is 
approximated by pUiUj and density fluctuations neglected in the acoustic source. 


Characteristics of the simulated turbulence 

The acoustic radiation from isotropic turbulence has been computed with the hybrid 
DNS method for the case with the following initial parameters 

Re x = 50, , Mt = 0.05 

where Re\ is the Taylor microscale Reynolds number, and Mt the turbulent Mach num- 
ber^ Note that Re\ = q\jv where q = \J u'u', A = q/\J w,- is the fluctuating 
vorticity, and u is the kinematic viscosity; while Mt = qjc where c is the mean speed 
of sound.) The initial energy spectrum is given by Eq. (3) with k m = 6. The initial 
data for the case discussed here is chosen to be incompressible; the velocity is solenoidal 
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(x = 0), the density is constant, and the pressure is initialized with the usual Poisson 
equation applicable to incompressible flows. Twenty simulations on a 64 3 grid with a 
time step of At — 0.00275 {K/t) 0 were caried out. A single, higher-resolution simulation 
was performed on a 128 3 spatial grid with At = 0.001375(A'/e) o , other parameters re- 
maining equal, to check the accuracy of the coarse-grid solution. The simulations were 
conducted for one eddy turnover time ( eot/Ko = 1 ). The 64 3 simulations were used 
to obtain ensemble-averaged statistics from the instantaneous acoustic pressure, while 
the 128 3 simulation was used to obtain the acoustic power from the appropriate flow 
statistics. 

The turbulence statistics at a given time in a DNS case are obtained by averaging over 
the computational volume. For example the turbulent kinetic energy K =< > /2 

where < . > denotes a volume average. In the present simulations, the turbulence is 
unforced. Consequently, the turbulent kinetic energy decays in time. Figs. 1 and 2 show 
the evolution of the turbulent kinetic energy K and turbulent dissipation rate e, respec- 
tively. K decreases by about a factor of 5 by the end of the simulation. The evolution of 
K for e 0 t/ Ko > 0.2 is well represented by a power law decay ( t — to)~ n with the exponent 
n — 1.4, a value which is slightly larger than experimentally measured exponents which 
lie in the range 1.15 < n < 1.35. The turbulent dissipation rate e increases during the 
initial transient in Fig. 2 due to the generation of small-scale fluctuations, and finally 
decays in the absence of external forcing. The microscale Reynolds number Re\ whose 
evolution is shown in Fig. 3 decays by a factor of 3 at the end of the simulation. The 
simulation is terminated at a nondimensional time tot/Ko = 1, because Re\ at later 
times could become too small to represent realistically the spatio-temporal complexity 
of turbulence. The skewness of the velocity derivative Sk = (du/dx) 3 /[(du/dx) 2 ] 3 ^ 2 is 
a measure of the non-linear vortex stretching. According to Lesieur(1990), experiments 
on grid turbulence give Sk ~ —0.4, while simulations give Sk ~ —0.5. Fig. 4 shows that 
skewness factor asymptotes to a value of Sk = —0.46 in the DNS which is consistent 
with previous experiments and simulations. 
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Characteristics of the computed acoustic pressure 

The fluctuating acoustic pressure is computed from the DNS data by the temporal 
form of the Lighthill analogy, 

(5) 

where [.] denotes evaluation at retarded time, while the overbar denotes ensemble average. 
Thus, P/i(x, t) is computed from Eq. (2); it’s ensemble average calculated over the 20 DNS 
cases; and finally, the ensemble average is subtracted out to calculate p^(x, t). Fig. 5 
shows the evolution of the acoustic pressure fluctuation p' A normalized by the ambient 
pressure P 0 . The smallness of the acoustic pressure fluctuation which is about seven 
orders of magnitude smaller than the thermodynamic pressure would present a formidable 
challenge for the direct computation of sound by simulation of the compressible Navier- 
Stokes equations in the far-field. Since the turbulent source decays with time, the acoustic 
pressure also decays with time. We have tried to factor out the decay in turbulence 
statistics by dividing the random pressure signal p' A (t ) by an appropriate deterministic 
function g(t). The choice of g(t ) was guided by the analysis of Proudman (1952) who 
considered the generation of noise by isotropic turbulence and used statistical models of 
various two-point moments within the framework of the Lighthill analogy to obtain the 
following expression for the acoustic power Pa per unit mass 

P A = a(u 3 //)^ (6) 

Co 

where a is a constant related to the shape of the longitudinal velocity correlation /(r), u 
is the root mean square (rms) of one of the velocity components, and l is the longitudinal 
integral length scale of the velocity. For the Heisenberg form of the energy spectrum, 
Proudman obtains 

Pa = 38(u 3 //)^ (7) 

C 0 

and, with the assumption /(r) = e - ^/*) 2 /^ Proudman obtains 

Pa = 13 (u 3 /l)4 (8) 

Co 
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Eq. (8) implies that g(t) oc yeM ( 5 can be used to normalize the acoustic pressure p a 
in order to factor out the temporal decay. Fig. 6 shows the evolution of rescaled acoustic 
pressure p* A defined by 

. = MMR. 

A y/(‘L/4)Ml 

Here, L is the length of the computational region. It appears that the normalization 
does factor out the temporal decay of the acoustic pressure. Fig. 7 shows the acoustic 
pressure fluctuations at different observer points. The instantaneous pressure signal at a 
given observation point is non-periodic, chaotic, and has a range of frequencies. 

The rms acoustic pressure Pmis(^) observed at a distance x is obtained by ensemble- 
averaging over 200 samples of fluctuating acoustic pressure p ' A . These samples are ob- 
tained from the 20 cases in the DNS with each case having 10 observation points per 
sphere distributed on two concentric spheres around the center of the flow domain. The 
acoustic power emitted from the volume V of fluid is 

P v (t) = r 2 

PoCo 

and the acoustic power per unit mass of turbulent fluid is 


M) = 


Mi 

PoV ' 


(9) 


Fig. 8 shows the acoustic power in decibels from the DNS and the theoretical result. 
The dotted curve in Fig. 8 corresponds to a value of a = 13 in Proudman’s work giving 

P A = 8.7«(()^ 

= 0.5 e(t)M t (t) 5 (10) 


and is obtained from Eq. (8) by using the definition of M t = \J2 K /cq and our DNS 
result e ~ 1.5 u 3 /l. The computed power is smaller than that obtained from Eq. (10) by 
6 — 8 dB over the time of the simulation. The ratio varies between 0.15 and 

0.07 in Fig. 9, i.e., by a factor of 2 during the decay of turbulent kinetic energy, while 
the numerator and denominator vary by about a factor of 15 during that time. Thus, 
the DNS is consistent with Proudman’s result that Pa is proportional to eM ( 5 . However, 
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if the proportionality constant from the DNS is taken to be P^/(M t 5 ~ 0.10, it is smaller 
relative to Proudman’s result, Eq. (10), by a factor of 5. 

Witkowska and Juve (1993) have computed the sound generated from unforced isotropic 
turbulence by performing large eddy simulations (LES) on 16 3 spatial grids for determin- 
ing the turbulent flow and then using the Lighthill acoustic analogy for computing the 
far-field sound. Five turbulent Mach numbers in the range of 0.012 to 0.015 were con- 
sidered; the initial Re\ was about 400. In agreement with our DNS results, the LES 
study finds that, although the acoustic power is proportional to cA/ ( 5 , it is 12 dB less 
than Proudman’s result for the Heisenberg spectrum. According to Witkowska and Juve 
(1993), a reconsideration of Proudman’s analysis with the von Karman spectrum instead 
of the Heisenberg spectrum decreases the theoretical estimate so that the LES result is 
6dB lower than the revised theoretical result. 

The turbulence in the DNS has moderate to low Reynolds number since Rex decays 
from a value of 50 to 15 during the simulations. The LES of Witkowska and Juve (1993) 
has a turbulence Reynolds number which is a factor of 10 higher than the DNS, albeit with 
the approximations inherent in a subgrid scale model. Although, there is a difference in 
Reynolds number, the DNS and LES studies are consistent in giving an acoustic power 
which is a factor of 5 smaller than Proudman’s result Eq. (10). Proudman’s analysis 
predates the present computational studies by about four decades and, of necessity, had 
to employ various assumptions such as Gaussian statistics for the velocity and its time 
derivatives, neglect of retarded time during evaluation of the integral expression for the 
acoustic power, neglect of the decay of unforced isotropic turbulence, and a specific 
shape for the longitudinal space correlation /(r). Lilley (1993), in a reconsideration of 
Proudman’s analysis, finds that the coefficient a in Eq. (6) is likely to be between 3 and 
10. In contrast, our DNS gives a ~ 2.6. The analysis of Lilley (1993) also elegantly 
shows that the acoustic power is sensitive to the shape of the space-time covariance of 
T xx , the flatness factor of the velocity, and is proportional to the fourth power of the 
characteristic Strouhal number s = ujI/u where l and u> are the integral scales associated 
with the spatial and temporal autocorrelations of the velocity, respectively, and u is the 
rms of a velocity component. Using the velocity flatness factor of 3.0 in our DNS results, 
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Lilley (1993) finds a = 3.4 with s = 1 and a = 8.3 with s = 1.25 in an analytical 
calculation which assumes the same /(r//) that corresponds to Proudman’s a = 13. 
Given the sensitvity of a to the assumed model correlations, it is not surprising that 
the analytically determined values for a are somewhat different from the corresponding 
results in the simulations. Overall, the DNS value of a is consistent with the lower end 
of the range of values of a found by Lilley (1993). 

The scaled acoustic power fU/eM ( 5 is significantly smaller than that typical of sub- 
sonic jets. In subsonic jets the acoustic efficiency (sound power /jet power) according to 
Lighthill (1954) is approximately 10~ 4 A4 5 , where M is the jet Mach number. The tur- 
bulence intensity u rms /Uj e t can be roughly estimated to be 0.17, and e is roughly 0.1 of 
the jet power, which implies that P^/eM, 5 ~ 7. Thus the sound radiated from isotropic, 
homogeneous turbulence is significantly smaller than that typical of jet acoustics. We 
argue below that one of the reasons for this dissimilarity is the homogeneity constraint. 

The acoustic pressure fluctuation is given by Eq.(5). Consider the integral in Eq. (5) 
as M t —* 0. 

j[Ta\-W]dy 
-* J Tij - Tij dy 

= V < >-V< % >. (11) 

The second line in Eq. (11) follows because retarded time can be neglected in the limit of 
M t —* 0, and the third line follows from the definition of the volume average denoted by < 
. >. In the case of homogeneous turbulence, < Tij > is equivalent to an ensemble average, 
and since further time averaging does not change a quantity which is already an average 

< > =< Tij >. Consequently the integral in Eq. (5) is zero when retarded time effects 
are neglected in the limit of M t — ► 0. However, in the case of inhomogeneous turbulence 

< Tij > is not a temporally smooth quantity, and this integral is not necessarily zero 
even if retarded time effects were neglected. 

For small but finite Mach number,, when retarded time effects cannot be neglected, 
it is clear that the volume integration smoothes out the temporal oscillation in [T,j] to a 
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greater extent in homogeneous turbulence than in inhomogeneous turbulence, which thus 
causes smaller fluctuations in the radiated sound from homogeneous turbulence relative 
to the inhomogeneous case. 

The frequency content of the acoustic pressure was obtained by analyzing the nor- 
malized acoustic pressure p A (t ) in Fourier space. The power spectrum E p (u > ) obtained by 
averaging over 400 records of p A (t) is shown in Fig. 10. The frequency is normalized by 
l 0 /u 0 , where / 0 is the initial longitudinal integral length scale and u 0 is the initial rms of 
a velocity component. The power spectrum in Fig. 10 has a peak at a Strouhal number 
St = ulo/uo — 3.5. Thus, the dominant frequency of the sound radiated to the far- field 
in the isotropic turbulence DNS is somewhat larger than the dominant frequency uq/Iq 
of the energy-containing range of turbulence. 

Computation of the acoustic power from the turbulence statistics 

The acoustic power obtained in our simulations is somewhat smaller than predicted 
by Proudman’s analysis. This discrepancy motivated an alternative method to obtain 
the acoustic power wherein the the fourth-order, two-point correlation T^ x {0)T^ s (r) from 
the DNS is used. 

Consider an observer point at location (x,0,0). Then, due to Eq. (5) for the acoustic 
pressure fluctuation, the acoustic power P A per unit mass at the observer point given by 
Eq. (9) becomes 

P A {t) = J dx A dx B (12) 

where V is the source volume, and and t , 4 denote the position of source point A and 
time of sound emission at point A so as to reach observer at time t\ similarly, x# and ts 
refer to source point B. Proudman assumes that T^. x {xA^t)T^ x (yiByt) is constant over the 
maximum time difference 50 that retarded time effects can be ignored during the 

volume integration with respect to source point B. This assumption leads to 

Pa ( t) = ^5 / TUxMM) dx A . (13) 
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Proudman, then assumes that t a can be considered constant during the integration over 
source point A, to obtain finally 


PA(t)=-^Jfl,(XA)fi,(x B )dxA (14) 

The two-point spatial autocorrelation in Eq. (14) is a function of r = |x^— x#| only, due to 
the turbulence isotropy. Let us define a\ = ( u 2 — < u 2 >,v 2 — < v 2 >,w 2 — < w 2 >), and 
denote the spatial correlation ( d 2 a' i /dt‘ 2 ) A (d 2 a , j /dt 2 ) B by <2,j(r). The subscript i = 1,2, 3 
denotes the x,y,z components, respectively. Then, Eq. (14) becomes 

p * m = i hJ Q " {r)dr (15) 

The following expression for Qij can be written due to the isotropy constraint (e.g. see 
Hinze(1975)) 

Qij(r) = (dV/diY [ f(r) ~ G(r) r i r j + G(r)«y] 

where the fluctuation (d 2 u 2 / dt 2 ) 1 = d 2 u 2 / dt 2 — < d 2 u 2 jdi 2 >, while F(r ) and G{r ) are 
the longitudinal and lateral correlation functions associated with Qij respectively. Thus, 

<2ii(r) = (d 2 u 2 /dt 2 )' 2 [F(r) cos 2 (6») + G(r) sin 2 (0)] 


where 6 is the angle between r and the x axis. The volume integration in Eq. (15) with 
respect to the separation vector r is then performed using a spherical coordinate system 
to yield 


Pa{*) = 


{d 2 u 2 /dt 2 ) 


2 V 2 rR r 2 


1 r j(F(r) + 2Gf(r))* 


(16) 


Proudman (1952) uses the statistical theory of isotropic turbulence to model the corre- 
lations in Eq. (16); we, on the other hand, directly compute ( d 2 u 2 /dt 2 )' 2 , F(r) and G(r ) 
from the DNS data base. 

The longitudinal correlation functions F(r) associated with Qij , and /(r) associated 
with the velocity field are shown in Fig. 11 at the time e 0 t/K 0 = 0.7. Evidently, the 
presence of the second time derivative leads to a much smaller correlation length of 
the acoustic source relative to the velocity. To confirm our results, a new simulation 
was performed with twice the resolution in both space and time. Fig. 12 shows that the 
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higher resolution 128 3 simulation gives a result for F(r) which is similar to the coarse-grid 
simulation. The correlation length scale for the acoustic source is 0.10, while the integral, 
Taylor and Kolmogorov length scales are 0.58, 0.12, and 0.03, respectively. It thus appears 
that the correlation length of d 2 u 2 jdt 2 in our simulation of isotropic turbulence is smaller 
than the turbulence integral length scale and larger than the dissipative length scale. 

Eq. (16) was used to provide an alternative estimate of the acoustic power from the 
DNS. Fig. 13 shows that the result from Eq. (16) is consistent with the previous result 
obtained from statistical post processing of the instantaneous acoustic pressure. Since 
retarded times are neglected in simplifying the exact expression for the acoustic power 
Eq.(12) to obtain Eq. (16), it is not surprising that there is some difference between the 
two methods. 

Conclusions 

We have applied the hybrid DNS approach to the problem of sound generated by 
isotropic turbulence. The combination of a spectral DNS of isotropic turbulence and 
the Lighthill acoustic analogy is demonstrated to be a feasible approach for obtaining 
the instantaneous acoustic pressure. Although the time derivative form of the Lighthill 
acoustic analogy requires more computer memory for storage than the spatial derivative 
form, it is preferable due to a less stringent time step constraint. 

The numerically computed acoustic efficiency of isotropic turbulence is much smaller 
than that in subsonic jet experiments. The frequency spectrum of the acoustic pressure 
indicates that the dominant acoustic frequency is somewhat larger than the frequency 
of the energy-containing eddies as found in Proudman (1952) and Lilley (1993). The 
computed acoustic power agrees with the theoretically derived proportionality to eM t 5 
deduced by Proudman (1952); however the constant of proportionality is smaller than 
the analytical result. The acoustic power was computed from the DNS both by ensemble- 
averaging the instantaneous acoustic pressure, and also from the two-point, fourth-order 
correlation f xx (Q)T xx (x). The two methods for obtaining the acoustic power yield con- 
sistent results. 


14 



Proudman (1952) analytically obtained the following relation for the acoustic power, 

P A = a(u 3 /l)^ 

where a — 13 with the particular choice /(r) = for the longitudinal velocity 

correlation. In contrast, a ~ 2.6 in the present DNS. The analysis of Lilley (1993) 
indicates that a is likely to have a value between 3 and 10. The analytical studies have 
had to necessarily make some assumptions about the space-time statistics of the acoustic 
source while the DNS, although free from such assumptions, is restricted to a moderate 
to small turbulence Reynolds number. The LES study of Witkowska and Juve (1993) 
which have turbulence Reynolds numbers larger than those in the DNS by an order of 
magnitude give a ~ 2.5. Although, the DNS and LES studies give consistent results, 
further DNS studies at higher Reynolds number are desirable in the future. 
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Fig. 3. Evolution of the microscale Reynolds number Re* in the DNS. 



Fig. 4. Evolution of the velocity derivative skewness in the DNS. 
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Fig. 7. Time series of the normalized acoustic pressure 
at four observation points in the DNS. 
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Fig. 8. The average power of the radiated sound. 
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Fig. 9. The scaled acoustic power in the DNS. 








Fig. 11. DNS results at time € Q t/K 0 =0.7 
on longitudinal spatial corralationa. 
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Fig. 12. Effect of increased spatial and temporal resolution on F(r), 
the longitudinal spatial correlation of the acoustic source. 
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Fig. 13. Comparison of two methods of obtaining the scaled acoustic power. 
Solid line is an ensemble average from the computed acoustic pressure, 
and symbols are from the computed two— point correlation of Tq . 
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